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Abstract 

Paddlefish  can  be  distinguished  from  other  freshwater  fish  by  the  presence 
of  a  long  paddle  shaped  snout  (Rostrum).  The  Rostrum  has  desirable 
mechanical  characteristics  that  could  be  explored  to  design  bio-inspired 
materials  that  have  better  strength  to  weight  ratio.  The  Rostrum  is  large 
and  generates  lift  like  an  aircraft  wing.  It  acts  as  a  stabilizer  preventing 
drag  during  filter  feeding  when  the  mouth  is  wide  open.  This  report  will 
discuss  the  Fluid  Structure  Interaction  simulations  conducted  to  analyze 
the  lift  generated  by  the  Rostrum,  fluid  flow  around  the  fish,  pressures  on 
the  Rostrum,  fluid  induced  swimming  enhancements,  and  stabilizing 
properties  of  the  Rostrum.  The  report  analyses  and  discusses  the 
experimentally  measured  forward  body  velocity  of  Paddlefish  during  filter 
feeding.  Laminar  flow  with  different  angles  of  attack  (o,  5,  and  10)  along 
and  against  the  direction  of  the  longitudinal  axis  of  the  swimming  fish  are 
also  studied  and  discussed. 
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1  Introduction 

1.1  Bio-inspiration 

Nature  has  evolved  from  a  single  cell  to  complex  biomechanics.  Several 
novel  materials  with  superior  properties  already  exist  in  nature.  For 
example,  the  woodpecker’s  head  has  been  a  source  of  inspiration  for 
designing  shock  absorbers.  Silk  from  silkworms  and  orb-weaving  spiders 
have  been  studied  to  incorporate  their  desirable  mechanical  properties  for 
use  in  biomaterials  (Vepari  and  Kaplan  2007).  One  property  that 
distinguishes  the  naturally  occurring  materials  from  man-made  materials 
is  the  unique  hierarchical  structure.  It  is  the  presence  of  this  hierarchical 
structure  that  imparts  these  materials  their  superior  strength  and 
performance  capability.  In  comparison  to  man-made  materials,  naturally 
occurring  materials  have  a  very  limited  number  of  basic  components. 
Advancement  towards  hierarchical  structure  gives  them  their  uniqueness, 
which  enables  them  to  withstand  adversity,  adapt,  and  heal.  Looking  at 
nature  from  an  engineer’s  eye  and  reverse  engineering  the  reason  for  the 
development  of  the  unique  hierarchical  structure  may  lead  us  towards 
designing  materials  that  have  broader  use  limits  and  areas  (Sen  2011). 

1.2  Paddlefish  Characteristics 

The  Rostrum  of  Paddlefish  is  a  unique  structure  comprising  of  a  network  of 
cartilage,  tissue,  and  interlocking  star  shaped  bones  called  stellate  bones. 
Figure  1  shows  the  stellate  bone  arrangement  in  the  Rostrum  of  Paddlefish. 
The  Rostrum  imparts  sensory,  electro  and  tactile  (Hoover  2013),  hydro¬ 
dynamics  and  strength  properties  depending  on  the  age  of  the  fish. 
Paddlefish  are  known  to  feed  efficiently  in  both  laminar  and  turbulent  flows. 
The  sensory  functions  of  the  Rostrum  enable  them  to  detect  the  type  of  flow 
(Gurgens  at  al.  2000).  Also,  the  sensory  functions  allow  them  to  detect  tiny 
zooplankton  by  detecting  them  void  of  sight  (Wilkens  et  al.  1997).  Further¬ 
more,  Paddlefish  readily  try  to  feed  on  dipoles  that  emitt  electrical  signals 
(Wilkens  and  Hofmann  2007).  Paddlefish  swim  with  their  mouth  wide  open 
during  filter  feeding  and  take  in  enormous  amounts  of  water.  The  gill  rakers 
help  them  filter  the  tiny  zooplankton  from  the  water.  During  filter  feeding, 
Paddlefish  swim  faster.  The  forward  body  velocity  is  used  to  transport  the 
water  at  high  speed  and  filter  food  quickly  and  efficiently.  The  function  of 
the  Rostrum  changes  at  different  stages  of  the  life  of  Paddlefish. 
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Figure  1.  Stellate  bone  of  Rostrum  (Hoover  2013). 


In  juvenile  stage,  the  shape  of  the  Rostrum  is  linear  and  size  and  area  of  the 
Rostrum  is  almost  one-third  the  body  length.  The  primary  function  of  the 
Rostrum  at  juvenile  stage  is  sensory.  During  the  sub-adult  stage  the  shape 
of  the  Rostrum  is  spatulate  and  the  primary  function  of  the  Rostrum  is 
hydrodynamic.  During  this  stage,  the  Paddlefish  are  active  filter  feeders.  In 
the  adult  stage,  the  shape  of  the  Rostrum  is  linear  and  the  primary  function 
is  mechanical.  Early  researchers  thought  the  Rostrum  was  used  to  dig  for 
food  (Forbes  1878, 1888a,  and  1888b;  Jordan  and  Evermann  1896; 
Stockard  1907;  Alexander  1914)  or  dislodge  organisms  from  vegetation 
(Norris  1923;  Beach  1902).  This  notion  was  put  aside  when  Paddlefish  were 
found  to  be  filter  feeders  (Sprague  1929;  Meyer  i960).  Adult  Paddlefish  can 
survive  without  Rostrum,  but  they  feed  less  efficiently  and  are  compara¬ 
tively  thinner.  Figure  2  shows  a  Paddlefish  with  Rostrum  injuries  (Hoover 
2013)  Paddlefish  are  extremely  mobile  species.  They  have  been  known  to 
travel  more  than  2000  miles  in  river  systems.  The  size,  shape,  and  position 
of  the  Rostrum  help  generate  lift  as  in  a  journal  article  written  by  Allen  and 
Riveros  (2013),  currently  in  review.  During  filter  feeding,  the  added  weight 
of  the  enormous  amount  of  water  would  pose  stability  issues  in  absence  of 
the  lift  generated  by  the  Rostrum. 
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Figure  2.  Rostrum  Injuries  (Hoover  2013). 


1.3  Paddlefish  population  decline 

Paddlefish  are  among  the  most  endangered  groups  of  fishes  in  North  America. 
They  were  once  abundant  in  most  central  US  river  systems.  Figure  3  shows 
the  population  data  of  Paddlefish  (Jennings  and  Zigler  2009).  One  of  the 
main  reasons  for  their  population  decline  is  overharvesting  (Williamson 
2003).  They  are  commercially  exploited  for  their  eggs,  or  roe  (caviar)  (Tower 
1908;  Stockard  1908;  and  Saffron  2002).  Additionally,  sedimentation  and  river 
modifications  are  other  reasons  for  their  population  decline.  One  of  the 
objectives  of  this  research  is  to  create  awareness  about  this  unique  fish  and 
learn  from  the  structural  and  mechanical  behavior  of  the  Rostrum  to  design 
novel  bio-inspired  material  that  would  have  better  strength  to  weight  ratio. 
These  bio-inspired  materials  could  have  potential  uses  in  the  area  of  body 
and  vehicle  armor,  ship  designs,  blast  loading  energy  dissipation,  repairs  of 
fatigue  and  fracture  in  navigation  steel  structures,  among  others. 

1.4  Objective 

Understanding  the  dynamic  effects,  forces,  and  pressures  involved  in  a 
swimming  fish  will  enhance  the  knowledge  in  various  applied  engineering 
fields,  i.e.  designing  submarines,  ships,  body  armors  etc.  Bergmann  and 
Iollo  (2011)  have  studied  the  flow  and  vortex  shedding  in  2D  and  3D  fish 
like  objects  that  resemble  an  airfoil.  Their  study  focuses  on  the  power 
required  by  the  fish  for  swimming,  both  for  a  single  fish  and  for  small  fish 
schools.  The  present  study  follows  the  experimental  work  done  by 
Sanderson  et  al.  (1994)  where  the  buccal  flow  velocity  of  the  Paddlefish  is 
measured  during  ram  ventilation,  ram  suspension  feeding,  and  prey 
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processing.  This  work  has  measured  the  forward  body  velocities  of 
Paddlefish  kept  in  stagnant  water  so  that  the  velocity  they  measure  is  the 
body  velocity  of  fish  and  not  the  water.  The  simulations  carried  out  in  the 
present  study  are  done  at  the  Reynolds  number  30,  proportional  to 
velocities  in  laminar  flow  regime  and  have  been  used  for  model  validation. 


Figure  3.  Paddlefish  population  distribution  (Jennings  and  Zigler  2009). 


■  STABLE 


NO  DATA 


Figure  4  shows  a  pictorial  representation  of  the  flowchart  depicting  the 
process  of  application  of  this  research.  It  shows  the  process  involved  in 
designing  novel  high  performance  materials.  The  end  product  of  such  cycle 
could  be  a  material  that  has  better  strength  to  weight  ratio  and  could  be 
used  to  make  lightweight,  high  performance  body  armor  for  our  soldiers. 

The  objective  of  this  study  can  be  classified  as  follows: 

1.  Conduct  fluid  structure  interaction  simulation, 

2.  Study  the  effect  of  flow  and  fish  velocity  on  the  Rostrum, 

3.  Observe  the  pressure  distribution  on  the  top  and  bottom  surface  of  the 
Rostrum,  and 

4.  Analyze  the  effect  of  angle  of  attack  on  the  uplift  pressure  and  flow 
velocity. 
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1.5  Fluid  Structure  Interaction  Simulations 

The  recent  advancements  of  computing  power  have  had  a  remarkable 
influence  on  the  nature  of  simulations  being  carried  out  on  high 
performance  computers.  Complicated  phenomenon  like  the  interaction  of 
one  or  more  solid  structure  with  the  surrounding  fluid  can  be  carried  out 
with  the  use  of  sophisticated  well-tested  algorithms.  The  effect  of  sub¬ 
merged  solid  structure  on  the  flow  pattern  can  be  studied  in  great  detail 
with  the  use  of  Multiphysics  techniques  like  the  fluid  structure  interaction 
(FSI).  Six  test  cases  are  simulated  for  the  present  study.  Laminar  flow  with 
different  angles  of  attack  (o,  5,  and  10)  and  velocities  of  the  swimming  fish 
are  studied  and  discussed.  The  report  discusses  the  lift  generated  by  the 
Rostrum,  fluid  flow  patterns  around  the  fish  body,  pressures  on  the 
Rostrum,  fluid  induced  swimming  enhancements,  and  stabilizing  properties 
of  the  Rostrum. 


Figure  4.  Bio-inspiration  flowchart. 
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2  Fluid  Structure  Interaction  Formulation  to 
Study  Flow  around  the  Polyondon 
Spathula  (Paddlefish) 

2.1  Multiphysics  Methods 

A  wide  range  of  interdisciplinary  problems  can  be  approached  using  fluid 
structure  interaction  (FSI).  With  the  recent  advancement  of  computer  power, 
a  variety  of  application  areas  can  be  addressed.  FSI  is  a  Multiphysics,  non¬ 
linear  problem  that  involves  the  interaction  between  incompressible  fluid 
flows  and  immersed  flexible  or  rigid  body  structures  (Hou  et  al.  2011).  FSI 
problems  analyze  the  interaction  of  one  or  more  solid  structures  with  the 
surrounding  fluid  flow.  FSI  problems  can  be  solved  with  two  approaches; 
monolithic  and  partitioned.  In  monolithic  approach,  the  fluid  and  structure 
dynamics  are  treated  in  a  single  mathematical  framework.  A  single  system 
equation  is  formed  for  the  whole  problem.  The  drawback  to  using  monolithic 
approach  is  the  amount  of  effort  that  goes  into  developing  specialized  code 
and  also  the  amount  of  labor  needed  to  maintain  such  specialized  code.  The 
partitioned  approach  treats  the  fluid  and  structure  model  separately  in 
different  mathematical  framework  as  the  name  suggests.  Existing  state  of  the 
art  solvers  are  used  to  address  the  fluid  and  structure  model.  Data  were 
exchanged  between  the  fluid  and  structure  solvers  at  the  interface  boundary. 
Although  this  approach  has  the  advantage  of  using  the  existing  legacy  codes, 
keeping  up  with  the  interface  boundary  which  changes  in  time  is  a  challenge. 
Partitioned  approach  is  used  in  the  present  study. 

The  present  study  uses  Abaqus  to  carry  out  the  FSI  simulations.  Various 
application-based  approaches  are  available  to  address  Multiphysics  prob¬ 
lems  (Table  l).  Multiphysics  approaches  with  abbreviations  are  listed  below: 

•  (CEL)  Coupled  Eulerian-Lagrangian  analysis, 

•  (SPH)  Smoothed  Particle  Hydrodynamics,  and 

•  (CFD/FSI)  Coupling  between  fluid  and  structure  model. 

FSI  methods  have  been  used  for  coupled  fluid  structure  interaction 
problems  in  biomechanics  applications  where  biofluids  such  as  blood 
interacts  with  biological  structures  such  as  arterial  walls  (Prasad  et  al. 
2011). 
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Table  1.  Application-based  comparison  of  Multiphysics  methods. 


CEL 

SPH 

CFD(FSI) 

Dam  break 

Blast  loading,  UNDEX 

V 

Impact  and  penetration  (Element  conversion) 

V 

Heat  Transfer 

Incompressible  fluid  flow  around  deformable  structures 

2.2  Computational  Fluid  Dynamics  (CFD)  Analysis 

The  integral  form  of  the  momentum  equation  for  an  arbitrary  control 
volume  can  be  written  as  the  following: 


— JpV  dV  +  JpV®VndS=-fVpdV  +  j'rndS  +  ffdV  (1) 


where: 


Vol  =  arbitrary  control  volume  with  surface  area  S, 
n  =  outward  normal  to  S, 
p  =  density  of  fluid, 

V  =  velocity  vector, 

V®V  =  tensor  product  of  velocity  vectors, 


u2  uv  uw 
V  <g>  V  =  vu  v2  vw 
wu  ivv  iv2 


(u,  v,  w) 
F 

T 


components  of  velocity  in  (x,  y,  z )  directions, 
body  force,  and 
viscous  shear  stress. 


Since  the  flow  is  incompressible, 


(2) 


V-V  =  0  (3) 

Rewriting  the  governing  equations  in  terms  of  the  primitive  or  physical 
variables  p,u,v,w,  andp,  we  have 
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‘du  du  du  du ^ 

- h  u - \-v - h  w — 

dt  dx  dy  dz 


dv  dv  dv  dv 

- h  U - \-V - \-lV - 

dt  dx  dy  dz 


dw  dw  dw  dw 
f  u - \-v - \-w 


dt  dx 


dy 


dz 


dp 


dp 

dy 

dp 

dz 


+  p 


+  9 


'  d2u 

d2u 

d2u 

dx2  ^ 

dy 2  + 

dz2 

/ 

d2v 

I 

d2v 

I 

d2v 

dx2 

\ 

dy2 

dz2 

/ 

d2w 

d2w 

_L 

d2 1 
_| _ 

dx2 

CM 

CO 

b 

dz 

+pgx 


+  P9n 


+pgz 


(4) 

(5) 

(6) 


Gravity  has  been  taken  into  consideration  as  a  body  force  (Toro  i999).The 
continuity  equation,  after  taking  the  incompressibility  into  account, 
reduces  to  the  following  form, 


du  ^  dv  ^  dw 
dx  dy  dz 


(7) 


We  have  four  equations  (4,  5,  6,  and  7)  and  four  unknowns  namely  u,v,w, 
and  p.  We  have  a  closed  system  of  equations  that  can  be  numerically 
solved  for  the  unknown  variables. 


2.2.1  Numerical  Implementation: 

Figure  5  shows  the  flowchart  of  the  incompressible  fluid  solver.  The  diver¬ 
gence  free  velocity  condition  imposed  because  of  incompressibility  makes 
the  numerical  solution  of  Navier-Stokes  Equation  a  challenge.  To  overcome 
this,  projection  method  is  used.  Projection  method  divides  the  velocity  field 
into  divergence  free  (solenoidal)  and  curl  free  (irrotational)  components 
using  Helmholtz  decomposition.  The  projection  algorithm  divides  the 
computation  of  velocity  in  two  stages.  In  the  first  stage,  an  intermediate 
velocity  that  does  not  satisfy  the  incompressibility  constraint  is  calculated. 
In  the  second  stage,  pressure  is  used  to  project  the  intermediate  velocity 
onto  a  space  of  a  divergence  free  velocity  field  to  obtain  the  updated  values 
of  pressure  and  velocity.  Hence, 


U=Usol+Uirrot=Usol+V<\> 


(8) 
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Figure  5.  Flowchart  of  the  incompressible  CFD  fluid  solver. 


Since,  V  x  Vcf>  =  0  for  some  scalar  functions  (<f>),  the  divergence  of  the 
equation  results  in  the  following: 


V-u  =  V2cf) 


Since, 


V-uso/=  0 


(9) 


(10) 


This  is  a  Poisson  equation  for  the  scalar  function^) .  If  the  vector  field  u  is 
known,  the  above  equation  can  be  solved  for  the  scalar  function(<f>)  and 
the  divergence  part  of  u  can  be  obtained  from  the  following  relation: 

usol=u-V  ^  (11) 


2.2.2  Least-square  gradient  approximation: 

The  discretization  of  viscous  terms  requires  gradient  of  primitive  variables 
across  the  cell  interfaces.  The  solution  method  uses  a  linearly  complete 
second  order  accurate  least-square  gradient  estimation.  This  facilitates 
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accurate  computation  of  dual  edge  fluxes  for  both  advective  and  diffusive 
processes. 

2.2.3  Advection  methods 

The  solver  uses  edge-based  treatment  for  advection  terms.  It  preserves  the 
smooth  variation  to  2nd  order  in  space.  The  advection  algorithm  uses  least- 
square  gradient  approximation  with  topologically  independent  slope 
limiters.  Sharp  gradient  are  captured  within  approximately  2-3  finite 
elements.  The  solver  uses  the  slope  limiters  in  addition  to  the  local 
diffusive  limiter  thereby  preventing  over/under  shoots  in  advected  fields 

2.3  Computation  Solid  Mechanics  Analysis 

The  explicit  dynamic  analysis  procedure  uses  explicit  integration  rule 
coupled  with  the  use  of  diagonal  “lumped”  element  mass  matrices.  Explicit 
central  difference  integration  rule  is  used  to  integrate  the  equation  of 
motion  of  the  body. 


£7<4)=£7'4>. 


A  t(i+1)  +Af(0 


U(i) 


(12) 


U(i+1)  =[/(i)  +A  t{i+1)U('+2)  (13) 

where  U  is  the  velocity  and  0  is  the  acceleration.  The  computational 
efficiency  associated  with  the  central  difference  procedure  is  because  of 
the  use  of  diagonal  element  mass  matrices.  The  inversion  of  the  mass 
matrix  used  for  computing  accelerations  is  triaxial. 

I7(0  /®)  (14) 


where  M  is  the  diagonal  lumped  mass  matrix,  F  is  the  applied  load  vector. 
And  I  is  the  internal  force  vector.  The  explicit  procedure  does  not  need 
iterations  and  tangent  stiffness  matrix. 

Special  attention  has  to  be  given  to  the  manner  in  which  mean  velocities  are 
handled  for  the  initial  condition.  For  presenting  results,  the  state  velocities 
are  stored  as  a  linear  interpolation  of  the  mean  velocities.  Figure  6  shows 
the  flowchart  of  the  numerical  implementation  of  the  solver  (Schutte  et  al. 
2010). 
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Figure  6.  Flowchart  of  explicit  structure  solver. 
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3  Co-Simulation  CFD/CSM  details 

3.1  Analysis  Details 

Several  methods  are  available  to  carry  out  Multiphysics  simulations  i.e. 
Coupled  Eulerian  Lagrangian  (CEL),  Smoothed  Particle  Hydrodynamics 
(SPH),  and  Fluid  Structure  Interaction  (FSI).  SPH  methods  are  best  suited 
for  applications  that  involve  extreme  deformation  (damage  models  and 
particle  movements)  like  a  blast  loading.  CEL  methods  are  best  suited  for 
applications  like  dam  break,  which  involves  the  interaction  between 
deformable  structures  with  compressible  multiphase  flows.  For  applications 
that  involve  the  interaction  of  fully  immersed  solid  deformable  body  with 
the  surrounding  incompressible  fluid,  FSI  methods  are  well  suited.  Based 
on  the  application  in  this  study,  FSI  method  was  selected.  Partitioned 
approach  is  used  for  the  fluid  structure  interaction  analysis.  The  equations 
governing  the  fluid  flow  and  displacement  of  structure  are  solved  separately. 
The  application  exchanges  data  (Boundary  conditions),  each  after  con¬ 
verged  iteration.  The  exchange  of  data  takes  place  at  the  fluid  structure 
interaction  boundary.  Rostrum  is  selected  for  exchanging  data  between  the 
two  solvers  since  it  is  a  deformable  body  and  the  structural  mechanical 
response  to  the  fluid  flow  is  of  interest  in  this  study.  Additionally,  there  is 
experimental  material  data  available  for  the  Rostrum  (Allison  at  al.  2013). 
The  analysis  considers  the  fish  body  as  a  display  body.  By  setting  the  fish 
body  as  a  display  body,  we  are  reducing  the  amount  of  calculations  done  on 
the  structure  solver  and  capturing  the  effects  of  the  body  on  the  flow 
surrounding  the  fish.  CFD  solver  considers  the  body  of  the  fish  for 
simulating  the  flow.  The  structure  model  does  computations  only  on  the 
Rostrum. 

Figure  7  shows  a  pictorial  representation  of  the  implementation  of  the 
native  FSI  capability  for  the  present  study.  Stabilized  segregated  approach  is 
used  in  which  the  structure  and  fluid  equations  are  solved  independently. 
The  interface  loads  and  boundary  conditions  are  exchanged  between  the 
solvers  after  a  converged  increment.  The  co-simulation  engine  handles  all 
the  interaction  of  loads  and  boundary  conditions  between  the  time 
dependent  fluid  and  structure  domain  without  any  user  intervention. 
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3.2 


Figure  7.  Partitioned  fluid-structure  interaction. 


CFD  Model  Details 


The  CFD  model  is  comprised  of  123506  linear  tetrahedral  elements  and 
24071  nodes.  The  total  number  of  variables  in  the  CFD  model  is  72213. 
Figure  8  shows  the  dimensions  of  the  CFD  model  used  for  the  simulations. 


Figure  8.  CFD  model  dimensions. 
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As  seen  in  Figure  8,  the  CFD  box  is  250  inches  long,  100  inches  deep,  and 
50  inches  wide.  The  geometry  of  the  fish  is  cut  out  from  the  geometry  of  the 
fluid  box  by  using  a  Boolean  operation  on  the  two  geometries  in  Abaqus. 
The  fish  is  placed  at  approximately  25  inches  in  depth  and  50  inches  away 
from  the  inflow.  The  boundary  conditions  are  discussed  in  chapter  4. 

3.3  CSM  Model  Details 

Material  properties  used  for  the  Rostrum  are  taken  from  the  experimental 
data  collected  by  Allison  et  al.  (2013)  and  Mississippi  State  University. 
Nano-indentation  was  carried  out  on  the  Rostrum  at  several  locations  as 
shown  in  Figure  9.  Modulus  was  calculated  using  the  force-displacement 
curves.  The  material  properties  used  for  the  simulations  are  taken  from  the 
values  shown  in  Table  2.  The  Rostrum  was  partitioned  at  approximate  loca¬ 
tions  where  the  data  were  collected.  Six  material  properties  are  used  for  the 
Rostrum  based  on  these  experimental  values.  Looking  at  Figures  9, 10,  and 
Table  2,  the  material  properties  for  the  middle-section  Ml,  is  higher  com¬ 
pared  to  M2.  Figure  10  shows  the  cross-section  of  the  Rostrum  from  the  fish 
mouth  to  its  tip.  Herein,  Ml  corresponds  to  slides  12  through  14  and  M2 
corresponds  to  slides  6  through  8.  Ml  has  more  compact  cartilage  arrange¬ 
ment  making  Ml  stiffer  than  M2.  Similarly,  M3  corresponds  to  slides  3 
through  5  in  Figure  10.  M3  is  stiffer  because  of  the  stellate  arrangement  of 
bones  as  seen  in  Figure  1. 


Figure  9.  Experimental  Nano-indentation  data  (Hoover  2013). 
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Table  2.  Average  result  for  Modulus  from  force-displacement  curves. 


Left  Side  (Gpa) 

Middle  Section  (Gpa) 

1 

12.48  ±  0.80 

9.22  ±  0.7 

2 

11.50  ±  0.50 

7.91  ±  0.34 

3 

12.16  ±  0.21 

13.67  ±  0.48 

The  geometry  of  Rostrum  is  created  from  the  point  cloud  data  of  the 
Rostrum  of  an  adult  Paddlefish.  3DMAX  was  used  to  extract  a  coarser  STL 
file  from  the  point  cloud  data.  Solidworks  was  used  to  extract  the  simplified 
geometry  in  IGES  format.  The  body  of  the  fish  was  constructed  using  the  2D 
drawing  shown  in  Figure  11.  Based  on  the  Rostrum  size,  the  dimensions  of 
the  fish  body  were  derived  (Hoover  2013).  Then,  Gridgen  was  used  to  create 
the  fish  body.  All  the  parts  were  imported  in  Abaqus  and  used  to  create  the 
whole  fish  body. 

The  structure  model  (Figure  11)  has  10802  linear  tetrahedral  elements  and 
2805  nodes.  The  total  number  of  variables  in  the  structure  model  is  8415. 
Computations  are  carried  out  only  for  the  Rostrum  part  (Figure  11).  The 
body  of  the  fish  is  set  up  as  Display  body.  The  Rostrum  and  fish  body  are 
attached  using  tie  constraints  and  velocity  boundary  conditions  are  applied. 
Three  points  are  selected  on  the  Rostrum  surface.  The  fish  body  follows 
these  points  as  the  Rostrum  moves.  Figure  11  shows  the  fish  model  used  in 
the  simulation.  Figure  12  shows  the  dimensions  of  the  Rostrum  used  for  the 
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structure  model.  Figure  13  shows  the  guidelines  used  (Hoover  2013)  to 
measure  the  dimensions  of  the  Paddlefish.  Figure  13,  along  with  the 
experimental  data  sheet,  provided  by  Hoover  (2013)  in  Figure  14,  which 
gives  guidelines  on  the  size  and  shape  of  Paddlefish  as  a  reference  to 
construct  the  body,  fins,  and  tail  of  the  Paddlefish  used  in  the  current 
simulations. 


Figure  11.  Fish  body  used  in  CSM. 


Figure  12.  Rostrum  dimensions. 
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Figure  13.  Rostrum  structure  model  dimensions  by  Hoover  (2013). 


Figure  14.  Rostrum  structure  model  dimensions  by  Hoover  (2013). 
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4  FSI  Results 

Time  dependent  fully  coupled  fluid  structure  interaction  simulations  are 
carried  out  in  the  present  study  to  analyze  the  hydrodynamic  and  structural 
properties  of  Rostrum  of  the  Paddlefish.  For  implementing  these  fully 
coupled  simulations,  an  incompressible  fluid  solver  is  used  to  computa¬ 
tionally  solve  the  laminar  flow  around  the  entire  fish.  The  fluid  solver  is 
coupled  with  the  structural  response  of  the  Rostrum  in  the  explicit  structure 
solver.  These  simulations  include  the  effect  of  shear  forces  on  the  Rostrum 
and  also  encompass  the  effect  of  deformation  of  Rostrum  on  the  fluid  flow. 
The  deformation  of  the  Rostrum  changes  the  boundary  condition  of  the 
fluid  flow.  The  interface  boundary  is  explicitly  defined  during  the  mesh- 
making  phase. 

4.1  Test  Case  Details 

Six  FSI  simulations  are  carried  out  in  the  present  study.  The  description  of 
the  flow  parameters  and  angle  of  attack  used  for  the  test  cases  is  shown  in 
Table  3.  Three  angles  of  attack  and  two  activities  of  swimming  fish  are 
analyzed.  Three  parametric  analyses  are  carried  out  with  fish  moving  in 
the  direction  of  fluid  flow  and  three  with  fish  moving  against  the  direction 
of  fluid  flow.  Very  low  velocities  are  used  for  the  investigation  to  validate 
the  model.  Further  research  in  the  area  will  involve  additional  parametric 
analyses  with  varying  velocities  of  water  as  well  as  the  fish  to  study  the 
hydrodynamic  and  structural  response  of  the  Rostrum  to  the  fluid  flow. 


Table  3.  Test  Case  Velocities. 


Activity 

0  Degree  (in/sec) 

5  Degree  (in/sec) 

10  Degree  (in/sec) 

Fish  swimming  with  the 
fluid  flow  (2.5  in/sec) 

Vx  =  2.75 

Vy  =  0 

Vx  =  2.739 

Vy  =  0.2397 

Vx  =  2.70822 

Vy  =  0.4775 

Fish  swimming  against 
fluid  flow  (2.5  in/sec) 

Vx  =  2.75 

Vy  =  0 

Vx  =  2.739 

Vy  =  0.2397 

Vx  =  2.70822 

Vy  =  0.4775 

4.2  Simulation  time  and  HPC  information 

All  simulations  were  carried  out  on  HPC-Diamond  (SGI  Altix  ICE  8200). 
Two  nodes  were  used  and  a  total  of  10  cores;  8  cores  for  CFD  and  2  cores 
for  CSM  solver.  It  took  an  average  of  l-hour  wall  time  to  simulate  0.1  sec 
real  time.  All  the  test  cases  were  executed  for  a  7-hour  wall  time.  One 
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extended  simulation  was  run  for  90  hours  to  ensure  that  the  test  cases 
were  run  long  enough  to  capture  the  effects  of  the  flow  around  the  fish  and 
to  ensure  that  the  flow  disturbance  induced  by  the  fish  geometry  reached  a 
stable  flow  pattern.  In  the  extended  run,  changes  around  the  body  of  the 
fish  after  0.5  sec  of  real  time  were  not  observed.  Figure  15  shows  the 
picture  of  the  SGI  Altix  ICE  8200  (Diamond)  located  at  the  ERDC  DSRC. 


Figure  15.  Diamond.  SGI  ALtix  Ice  -  172  TFLOPS. 


4.3  Velocity  profile  for  fish  swimming  against  the  flow  direction 

Figure  16  shows  the  velocity  profile  for  fish  moving  against  the  direction  of 
flow.  The  angle  of  attack  used  for  this  execution  is  zero  degrees.  The  fish  is 
moving  in  an  opposite  direction  of  the  fluid  flow.  When  the  forward  moving 
Rostrum  velocity  comes  in  contact  with  the  opposite  flow  direction  velocity 
vectors;  tiny  vortices  are  created  near  the  top  and  bottom  surface  of  the 
Rostrum.  The  opposite  direction  of  velocity  vectors  creates  these  vortices.  A 
high  velocity  region  is  seen  very  close  to  the  upper  and  lower  surface  of  the 
Rostrum  as  seen  in  Figure  16.  At  some  distance  from  the  Rostrum  the 
magnitude  of  velocity  vectors  reaches  zero.  The  zero  velocity  cloud  is  seen  at 
this  point  as  seen  in  Figure  16.  After  the  zero  velocity  cloud,  the  velocity 
increases  steadily  until  it  matches  the  flow  velocity.  The  velocity  distribution 
above  and  below  the  surface  of  Rostrum  looks  symmetric  for  angle  of  attack 
of  zero  degree  as  seen  in  Figure  16. 

Figure  17  depicts  the  velocity  profile  for  fish  moving  against  the  direction  of 
flow.  The  angle  of  attack  used  for  this  run  is  5  degrees.  A  slight  increase  in 
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velocity  is  seen  near  the  tip  of  the  Rostrum.  As  seen  in  Figure  17,  an  uneven 
velocity  distribution  is  caused  at  the  tip  of  the  Rostrum  when  the  angle  of 
attack  is  not  zero.  This  helps  propel  the  fish  in  upward  forward  direction. 


Figure  16.  Velocity  profile  for  fish  moving  against  the  flow  direction  at  an  angle  of  attack  of 

0  degrees. 


Figure  17.  Velocity  profile  for  fish  moving  against  the  flow  direction  at  an  angle  of  attack  of  5  degrees. 


Figure  18  demonstrates  the  velocity  profile  for  fish  moving  against  the 
direction  of  flow  for  an  angle  of  attack  of  10  degrees.  An  uneven  distribution 
of  velocity  is  observed  at  the  tip  of  the  Rostrum.  A  similar  observation,  like 
that  of  Figure  17,  is  seen  for  the  angle  of  attack  of  10  degrees.  The  increase  in 
velocity  at  the  Rostrum  tip  is  more  pronounced  as  expected. 
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Figure  18.  Velocity  profile  for  fish  moving  against  the  flow  direction  at  an  angle  of  attack  of  10  degrees. 


Figures  16  through  18  are  the  velocity  distribution  profile  at  the  last  time 
step  of  the  7-hour  simulations.  The  velocity  was  stabilized  after  0.5  seconds 
of  simulation  time.  The  velocity  distribution  at  the  tip  of  the  Rostrum 
becomes  uneven  with  the  increase  in  angle  of  attack. 

4.4  Velocity  profile  for  fish  swimming  along  the  flow  direction 

When  the  fish  is  moving  along  the  longitudinal  axis  of  the  flow  direction,  the 
flow  does  not  offer  much  resistance  to  the  movement  of  the  fish.  Figure  19 
shows  the  velocity  profile  for  zero  angle  of  attack.  Very  uniform  velocity 
distribution  is  observed  above  and  below  the  Rostrum.  The  velocity 
increases  as  a  function  of  distance  as  we  move  away  from  the  body  of  the 
fish  until  it  matches  the  flow  velocity  as  expected. 


Figure  19.  Velocity  profile  for  fish  moving  along  the  flow  direction  at  an  angle  of  attack  of  0  degrees. 
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Figure  20  depicts  the  velocity  distribution  when  the  angle  of  attack  of  the 
inflow  is  5  degrees.  A  slight  difference  in  velocity  is  seen  above  and  below 
the  surface  of  Rostrum.  Velocity  increases,  as  a  function  of  distance,  as  we 
move  away  from  the  body  of  fish. 


Figure  20.  Velocity  profile  for  fish  moving  along  the  flow  direction  at  an  angle  of  attack  of  5  degrees. 


Figure  21  shows  the  velocity  profile  when  the  angle  of  attack  is  10  degrees. 
An  uneven  distribution  of  velocity  is  observed  at  the  tip  of  Rostrum  as 
expected. 


Figure  21.  Velocity  profile  for  fish  moving  along  the  flow  direction  at  an  angle  of  attack  of  10  degrees. 


Figures  19  through  21  are  the  velocity  profile  at  the  last  time  step  of  the 
7-hour  simulations. 
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4.5  Pressure  profile  for  fish  swimming  against  the  flow  direction 

From  the  basic  governing  laws  of  physics,  everything  in  nature  has  a 
tendency  to  evolve  to  a  state  where  survival  is  possible  with  the  least 
possible  expenditure  of  energy.  Fish  are  good  examples  of  this  process.  The 
body  of  a  fish  is  streamlined  so  that  it  experiences  minimum  stresses  and 
also  dissipates  energy  efficiently.  In  the  present  study,  we  have  studied  the 
pressure  distribution  above  and  below  the  surface  of  Rostrum. 

Figure  22  shows  the  pressure  contours  when  the  fish  is  moving  against  the 
flow  direction.  The  angle  of  attack  is  zero  degrees.  A  high-pressure  contour 
of  the  order  of  magnitude  io-3  is  seen  near  the  tip  of  the  Rostrum.  This 
high-pressure  contour  is  a  result  of  the  reaction  force  between  the  fluid 
and  fish  moving  in  opposite  directions.  A  small  high-pressure  bubble  is 
seen  near  the  mouth  of  the  fish.  This  high-pressure  bubble  will  help  push 
the  fish  upwards  and  assist  in  swimming.  High  pressures  are  also  seen 
when  the  flow  comes  in  contact  with  the  fins  and  tails  as  seen  in  Figure  22. 


Figure  22.  Pressure  profile  for  fish  moving  against  the  flow  direction  at  an  angle  of  attack  of  0  degrees. 
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Figure  23  represents  the  pressure  contours  when  the  fish  is  moving  against 
the  direction  of  flow.  The  angle  of  attack  is  5  degrees.  The  angle  of  attack 
generates  high-pressure  on  the  bottom  of  the  Rostrum  of  the  order  of 
magnitude  lOA  Paddlefish  swim  with  their  mouth  wide  open  during  filter 
feeding.  They  take  in  enormous  amounts  of  water  during  this  process  and 
transport  this  water  out  at  a  very  high  velocity.  The  forward  body  velocity  is 
used  to  transport  this  water.  The  shape,  size,  and  position  of  Rostrum 
generate  lift,  which  assists  the  fish  during  this  process  and  also  prevents 
nose-diving. 
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Figure  23.  Pressure  profile  for  fish  moving  against  the  flow  direction  at  an  angle  of  attack  of  5  degrees. 
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Figure  24  depicts  the  pressure  contours  when  the  fish  is  moving  against  the 
direction  of  flow.  The  angle  of  attack  is  10  degrees.  Figure  24  shows  the 
negative  pressure  pocket  generated  near  the  mouth  of  the  fish.  This  negative 
pressure  helps  in  creating  the  vacuum  that  helps  the  fish  gulp  in  water 
quickly.  High-pressure  contours  are  seen  at  the  bottom  of  the  Rostrum. 


Figure  24.  Pressure  profile  for  fish  moving  against  the  flow  direction  at  an  angle  of  attack  of  10  degrees. 
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As  seen  in  Figures  22  through  24,  the  magnitude  of  the  high-pressure  is 
around  the  same  range  for  all  angles  of  attack.  However,  when  the  angle  of 
attack  is  zero  degrees,  high-pressure  contour  is  seen  at  the  tip  and  mouth 
of  the  fish.  The  bottom  middle  section  of  the  Rostrum  area  does  not 
experience  the  high-pressure.  When  the  angle  of  attack  is  5  degrees,  the 


ERDC/ITL  TR-13-4 


25 


high-pressure  contour  covers  the  entire  bottom  area  of  the  Rostrum  and  is 
more  pronounced  on  the  lower  body  area  of  the  fish.  Figure  24  shows  that 
when  the  angle  of  attack  is  10  degrees,  the  high-pressure  contour  covers 
wider  areas  on  the  bottom  surface  of  the  Rostrum.  Also,  the  mouth  of  the 
fish  has  a  small  negative  pressure  pocket. 

4.6  Pressure  profile  for  fish  swimming  along  the  flow  direction 

Pressure  contours  are  created  for  the  simulations  when  the  fish  is  moving 
in  the  direction  of  flow  for  the  three  angles  of  attack.  Since  the  fish  is 
moving  in  the  direction  of  the  flow,  very  low  resistance  is  offered  by  the 
flow  to  the  movement  of  the  fish.  Hence,  the  reaction  force  and  pressures 
on  the  Rostrum  are  negligible. 


Figure  25  shows  the  pressure  contours  when  the  fish  is  moving  in  the 
direction  of  the  flow  at  an  angle  of  attack  of  zero  degrees.  Since  the  fish  is 
moving  in  the  direction  of  the  fluid  flow,  the  fluid  pushes  the  fish  in  the 
forward  direction.  This  will  generate  high-pressure  near  the  tail  and  back 
of  fins  as  seen  in  Figure  25. 


Figure  25.  Pressure  profile  for  fish  moving  along  the  flow  direction  at  an  angle  of  attack  of 

0  degrees. 


PRESSURE 

(Avg:  75%) 

r  +2.01Se-03 

-  -l-l .000e-03 

-l-5.000e-04 

+8.14Se-10 

-5 .000e-04 
-1.000e-03 
-1.500e-03 
-2 .000e-03 
-2 ,500e-03 
-3.000e-03 
-3 ,500e-03 

-4.000e-03 

-4.500e-03 

-  -5 .000e-03 

L  -2 ,385e-02 

ODB:  CFD-Box_0-5.odb  Abaqus/CFD  6.12-3  Wed  Apr  17  08:56:23  Central  Daylight 


,  y  Step:  Step-1,  0  deg  incline 

Increment  3201:  Step  Time  =  0.8201 

Primary  Var:  PRESSURE 

Deformed  Var:  U  Deformation  Scale  Factor:  +  1.000e+00 


Figure  26  depicts  the  pressure  profile  when  the  fish  is  moving  in  the 
direction  of  the  flow  at  an  angle  of  attack  of  5  degrees.  A  high-pressure 
region  is  seen  below  the  Rostrum.  The  angle  of  attack  coupled  with  the 
shape,  size,  and  position  of  the  Rostrum  helps  create  this  high-pressure 
region. 
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Figure  26.  Pressure  profile  for  fish  moving  along  the  flow  direction  at  an  angle  of  attack  of  5  degrees. 
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Figure  27  represents  the  pressure  distribution  when  the  fish  is  moving  in 
the  direction  of  flow  at  an  angle  of  attack  of  10  degrees.  The  increase  in 
angle  of  attack  makes  the  high-pressure  region  below  the  Rostrum  surface 
more  pronounced.  High-pressure  regions  below  the  surface  of  Rostrum 
helps  the  fish  swim  efficiently  during  filter  feeding.  The  size,  shape,  and 
position  of  the  Rostrum  help  stabilize  the  fish. 


Figure  27.  Pressure  profile  for  fish  moving  along  the  flow  direction  at  an  angle  of  attack  of  10  degrees. 
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Figures  25  through  27  show  the  high-pressure  contours  are  the  same 
magnitude  (to-10)  for  the  three  angles  of  attack.  However,  with  the  increase 
in  angle  of  attack  they  spread  to  larger  areas  below  the  Rostrum  of  the 
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Paddlefish.  Also,  since  the  flow  is  pushing  the  fish  from  behind,  higher 
pressures  of  magnitude  (ich)  are  seen  at  the  tail  and  back  of  fins. 

4.7  Pressure  graph  near  mouth  of  fish  for  the  three  angles  of  attack 

A  reference  point  was  picked  near  the  mouth  of  the  fish  as  seen  in  Figure  29. 
Pressure  was  plotted  as  a  function  of  time  for  the  three  angles  of  attack.  For 
all  three  cases,  the  fish  was  moving  against  the  direction  of  flow.  Negative 
pressure  peaks  are  observed  when  the  fish  is  moving  at  an  angle  of  attack  of 
ten  degrees  (Figure  28).  These  negative  pressure  peaks  help  in  the  filter 
feeding  process.  Since  these  peaks  are  generated  near  the  mouth  of  the  fish, 
they  help  in  creating  vacuum  near  the  mouth  of  the  fish,  which  in  turn 
assists  in  the  filter  feeding  process  (Sanderson  et  al.  1994). 

4.8  Pressure  on  fins  of  the  Paddlefish 

Figure  30  shows  the  pressure  on  the  four  fins.  For  fins  1  and  2,  the 
pressure  is  positive  at  the  beginning  of  the  simulation  as  seen  in  Figure  31. 
For  fins  3  and  4,  the  pressure  is  negative  at  the  beginning  of  the  simulation 
as  seen  in  Figure  32.  The  overall  pressure  distribution  follows  the  same 
trend  for  all  the  fins.  This  is  the  test  case  where  the  fish  is  swimming 
against  the  flow  direction.  The  angle  of  attack  is  10  degrees.  Negative 
pressure  peaks  are  also  observed  as  seen  in  Figures  31  and  32. 


Figure  28.  Pressure  near  the  mouth  offish  for  the  three  angles  of  attack. 
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Figure  29.  Pressure  reference  point  near  the  mouth  of  fish  for  the  three  angles  of  attack. 
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Figure  30.  Reference  point  on  fins  1,  2, 3,  and  4. 


4.9  Pressure  graph  on  bottom/top  surface  of  Rostrum 

Pressure  is  plotted  as  a  function  of  true  distance  along  path  on  the  top  and 
bottom  surface  of  the  Rostrum  (Figure  33  and  Figure  34).  Pressure  is 
plotted  as  a  function  of  distance  for  the  case  in  which  the  fish  is  moving  in 
the  direction  of  flow  for  the  three  angles  of  attack  o,  5,  and  10  degrees.  As 
seen  in  Figures  35-37,  very  low  pressures  are  seen  on  the  top  and  bottom 
surface  of  the  Rostrum.  Also,  the  pressures  on  top  are  slightly  lower  as 
compared  to  the  bottom  of  the  Rostrum.  As  seen  in  Figure  35,  when  the 
angle  of  attack  is  zero  degrees,  the  pressure  on  top  and  bottom  surface  are 
almost  identical.  As  the  angle  of  attack  increases,  we  see  higher  pressure 
being  developed  on  the  bottom  surface  of  the  Rostrum.  These  higher 
pressures  provide  lift,  which  in  turn  propels  the  fish  in  upward  direction 
and  also  prevents  nose-diving  to  bottom  during  filter  feeding. 
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Figure  32.  Pressure  on  fins  3  and  4. 
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Figure  33.  Path  reference  on  Rostrum  top. 


Figure  34.  Path  reference  on  Rostrum  bottom. 
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Figure  35.  Pressure  as  a  function  of  true  distance  on  Rostrum  top  &  bottom  at  an  angle  of 

attack  of  10  degrees. 


Figure  36.  Pressure  as  a  function  of  true  distance  on  Rostrum  top  &  bottom  at  an  angle  of 

attack  of  5  degrees. 
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Figure  37.  Pressure  as  a  function  of  true  distance  on  Rostrum  top  &  bottom  at  an  angle  of  attack 

of  0  degrees. 


4.10  Velocity  graph  at  Rostrum  tip  (1) 

An  uneven  distribution  of  velocity  is  seen  at  the  top  and  bottom  surface  of 
the  Rostrum  with  the  increase  in  the  angle  of  attack.  This  effect  is  more 
evident  at  the  tip  of  the  Rostrum.  Velocity  is  plotted  as  a  function  of  time 
to  see  the  effect  of  the  angle  of  attack  on  the  velocities.  Figure  38  shows  the 
reference  point  that  was  picked  near  the  tip  of  the  Rostrum.  Velocity  is 
plotted  as  a  function  of  time  for  the  case  when  the  fish  is  moving  against 
the  flow  at  the  three  angles  of  attack.  As  seen  in  Figure  39,  with  the 
increase  in  angle  of  attack,  the  velocity  at  the  reference  point  increases. 

The  angle  of  attack  creates  high  pressure  on  the  bottom  of  the  Rostrum  as 
seen  from  the  pressure  contour  plots.  The  high-pressure  on  the  bottom  of 
Rostrum  pushes  the  fish  in  an  upward  direction.  This  effect  coupled  with 
the  body  velocity  of  the  fish  is  responsible  for  the  increase  in  the  velocity  at 
the  Rostrum  reference  point  shown  in  Figure  38. 

4.11  Velocity  graph  at  Rostrum  tip  (2) 

Another  point  was  picked  at  the  tip  of  the  Rostrum  to  analyze  the  increase 
in  velocity  with  the  increase  in  angle  of  attack.  The  reference  point  shown 
in  Figure  38  is  chosen  very  close  to  the  tip  of  the  Rostrum.  Since  there  is  a 
lot  of  dynamics  involved  at  the  Rostrum  tip,  the  velocity  distribution  is 


ERDC/ITL  TR-13-4 


33 


studied  at  another  point  as  well.  This  reference  point  is  chosen  farther 
away  from  the  tip  of  the  Rostrum.  Figure  40  shows  the  reference  point 
picked  at  the  tip  of  the  Rostrum.  As  seen  in  Figure  41,  the  velocity  at  the 
reference  point  increases  with  the  increase  in  angle  of  attack.  The  zero 
velocity  cloud  is  not  seen  around  this  point.  The  lift  generated  because  of 
the  Rostrum  and  the  angle  of  attack  is  responsible  for  the  increase  in 
velocity  around  this  point. 


Figure  38.  Velocity  reference  point  at  Rostrum  tip  (1) 


Figure  39.  Velocity  as  a  function  of  time  at  Rostrum  tip  (1) 


2 

1.9 

1.8 

1.7 

^1.6 

<U  1  5 

to 

' - 

jE  1 .4 

>1  2 
^  1.3 

u 

-2  1.2  H 

:> 

1.1 

1 


0.1  0.2  0.3  0.4  0.5 

Time  (Sec) 


0.6 


0.7 


0.8 


ERDC/ITL  TR-13-4 


34 


Figure  40.  Velocity  reference  point  at  Rostrum  tip  (2) 


Figure  41.  Velocity  as  a  function  of  time  at  Rostrum  tip  (2) 


4.12  Streamlines  depicting  flow  patterns 

Streamlines  help  understand  and  visualize  the  flow  patterns  that  are  not 
visible  to  eyes  while  observing  a  physical  phenomena  occurring  in  nature. 
When  the  fish  is  swimming  against  the  direction  of  flow,  the  water  particles 
hit  the  body  of  the  fish  and  get  reflected  in  multiple  directions  creating 
vortices  above  and  below  the  Rostrum.  Figure  42  shows  the  flow  pattern 
when  the  fish  is  moving  against  the  flow  direction.  The  flow  shows  a  signifi¬ 
cant  amount  of  disturbance  above  and  below  the  surface  of  the  Rostrum.  As 
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the  flow  nears  the  body  of  the  fish,  the  vortices  and  turbulences  generated 
by  the  paddle  begin  to  dissipate,  developing  a  streamlined  flow.  Figure  42 
also  shows  how  the  Rostrum  helps  stabilize  the  flow  and  withstands  the 
disturbance  generated  when  the  fish  is  moving  against  the  direction  of  flow. 


Figure  42.  Streamlines  when  fish  is  moving  against  the  flow  direction 
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Figure  43  shows  the  flow  pattern  when  the  fish  is  moving  in  the  direction 
of  flow.  The  flow  is  very  streamlined  and  conforms  to  the  body  of  the  fish 
as  visible  in  Figure  43.  No  disturbances  were  seen  near  the  Rostrum.  The 
reaction  forces  and  pressures  experienced  by  the  Rostrum  in  this  case  are 
very  negligible  as  discussed  in  section  4.6  earlier.  In  this  case,  the  flow  of 
the  water  helps  the  fish  to  propel  in  forward  direction,  increasing  the 
velocity  of  the  fish  thereby  reducing  the  energy  and  muscle  movement 
required  to  swim 


Figure  43.  Streamlines  when  fish  is  moving  in  the  flow  direction 
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5  Conclusions  and  Recommendations 

5.1  Conclusions 

The  size,  shape,  and  position  of  the  Rostrum  develop  the  necessary  lift 
required  by  the  fish  for  filter  feeding  and  balance.  During  filter  feeding, 
Paddlefish  take  in  enormous  amounts  of  water.  This  added  weight  of  water 
would  make  it  extremely  difficult  for  the  Paddlefish  to  swim  efficiently. 
Also,  experimental  data  collected  on  filter  feeding  Paddlefish  show  that 
they  increase  their  body  velocity  by  ~6o  percent  to  transport  this  water  at 
very  high  speeds  and  filter  the  tiny  zooplanktons  from  water.  To  make  this 
process  more  efficient  with  minimal  amount  of  muscle  activity  from  the 
fish,  the  lift  generated  by  the  Rostrum  is  extremely  helpful  and  required. 
The  lift  generated  by  the  presence  of  Rostrum  will  reduce  the  amount  of 
muscle  movement  required  to  increase  the  forward  body  velocity.  Since 
the  shape  of  the  fish  is  streamlined,  the  Rostrum  as  well  as  the  fish  body 
experience  very  minimum  pressures.  With  the  increase  in  the  angle  of 
attack,  the  pressures  on  the  bottom  of  the  Rostrum  increase  thereby 
facilitating  the  forward  movement  of  the  fish  just  like  an  aircraft  wing. 

5.2  Velocity  of  fish 

Based  on  experimental  data  collected  for  Paddlefish  kept  in  round  tanks  in 
the  laboratory,  Paddlefish  increase  their  body  velocity  by  about  6o  percent 
during  filter  feeding.  Fish  move  by  contracting  muscles  that  propel  their 
body  in  the  forward  direction.  During  filter  feeding,  Paddlefish  take  in 
enormous  amounts  of  water  in  their  mouth.  Swimming  with  the  added 
weight  of  this  water  will  require  more  effort  from  the  Paddlefish.  In 
addition  to  this  work,  there  will  be  the  fact  that  Paddlefish  use  forward 
body  velocity  to  transport  this  water  at  high  velocity.  The  lift  generated  by 
the  Rostrum  during  this  phase  will  tremendously  reduce  the  amount  of 
work  required  by  the  fish  to  move  at  high  speed.  Also,  when  the  fish  is 
moving  against  the  flow  of  water,  Rostrum  acts  as  a  stabilizer.  It  forms  a 
zero  velocity  cloud  on  top  and  bottom  surface  of  the  Rostrum  that  takes 
care  of  the  vortices  generated  due  to  the  fish  body  and  makes  the  fluid  flow 
streamlined  when  it  reaches  the  body  of  the  fish. 


ERDC/ITL  TR-13-4 


37 


5.3  Pressure  on  Rostrum 

The  body  of  fish  is  streamlined  and  designed  to  withstand  the  static 
pressure  and  the  pressures  generated  due  to  flow  around  the  fish  body. 
Uplift  pressures  on  the  Rostrum  are  generated  to  aid  the  fish  swim 
efficiently  when  the  angle  of  attack  increases.  In  the  current  study,  we  have 
found  very  low  pressures  on  the  surface  of  the  Rostrum.  Also,  high  uplift 
pressure  is  generated  on  the  bottom  of  the  Rostrum  that  helps  generate 
lift.  Negative  pressure  pockets  were  seen  at  the  angle  of  attack  of  10 
degrees.  These  negative  pressure  pockets  create  the  vacuum  effect  that 
aids  the  fish  in  taking  in  large  amounts  of  water  during  filter  feeding.  The 
mechanical  property  and  shape  of  the  Rostrum  is  designed  such  that  it 
dissipates  energy  efficiently,  it  is  lightweight,  withstands  pressure  load 
efficiently,  and  enhances  the  generation  of  vortices  that  develop  uplift 
pressures  during  filter  feeding. 

5.4  Angle  of  attack  of  Rostrum 

The  increase  in  angle  of  attack  of  the  flow  generates  higher  pressure  on  the 
bottom  of  the  Rostrum.  These  increases  in  uplift  pressures  help  generate 
lift  that  propels  the  fish  in  the  forward  direction.  During  filter  feeding,  the 
fish  takes  in  enormous  amount  of  water  and  transports  this  water  at  very 
high  velocity.  Paddlefish  uses  forward  body  velocity  to  transport  this 
water.  The  lift  generated  by  the  Rostrum  in  conjunction  with  the  angle  of 
attack  helps  in  this  process.  We  have  observed  that  the  increase  in  angle  of 
attack  generates  negative  pressure  peaks  that  aid  in  filter  feeding. 

5.5  Recommendations  for  Further  Research 

1.  Include  the  hierarchical  structure  of  the  Rostrum  in  the  structure  model. 
Although  the  experimentally  measured  mechanical  properties  for  Rostrum 
are  used  for  the  structural  model,  the  incorporation  and  the  development 
of  hierarchical  structure  will  throw  light  towards  developing  better 
performing  materials. 

2.  Compare  different  Reynolds  number  flows.  Since  Paddlefish  are  known  to 
feed  in  both  laminar  and  turbulent  flows,  it  will  be  beneficial  to  see  the 
flow  pattern  around  the  fish  and  the  structural  properties  of  Rostrum  in 
turbulent  flows. 

3.  Since  the  main  focus  of  this  study  is  on  the  lift  generated  by  the  Rostrum 
during  filter  feeding,  including  the  mouth  opening  phenomena  during  ram 
filter  feeding  seems  to  be  the  next  logical  step. 
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4.  The  experimental  data  available  for  Paddlefish  during  filter  feeding  is  very 
limited.  Future  research  needs  to  concentrate  on  performing  experiments 
to  determine  fish  swimming  patterns  and  velocities,  pressures,  angles  of 
attack,  Reynolds  number  of  flow,  etc. 

5.  The  lack  of  material  data  on  fish  body  has  led  us  to  include  the  body  of  fish 
as  a  display  body.  Future  research  efforts  maybe  directed  towards 
including  the  muscle  movement  of  the  fish  body  in  the  structure  model  to 
simulate  a  more  realistic  swimming  fish. 
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